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ABSTRACT 

Magnetic field estimates for nearby isolated neutron stars (INS) help to constrain both the character- 
istics of the population and the nature of their peculiar X-ray spectra. From a series of XMM-Newton 
observations of RX J2143. 0+0654, we measure a spin-down rate of v = (—4.6 ± 2.0) x 10 _16 Hzs _1 . 
While this does not allow a definitive measurement of the dipole magnetic field strength, fields of 
> 10 14 G such as those inferred from the presence of a spectral absorption feature at 0.75 keV are 
excluded. Instead, the field is most likely around 2 x 10 13 G, very similar to those of other INS. We 
not only suggest that this similarity most likely reflects the influence of magnetic field decay on this 
population, but also discuss a more speculative possibility that it results from peculiar conditions on 
the neutron-star surface. We find no evidence for spectral variability above the ~ 2% level. We con- 
firm the presence of the 0.75-keV feature found earlier, and find tentative evidence for an additional 
absorption feature at 0.4 keV. 

Subject headings: stars: individual (RX J2143. 0+0654) — stars: neutron — X-rays: stars 



1. INTRODUCTION 



Th e so -called isolate d neutron stars (INS; see iHaberj 
2007] and iKaplanl 120081 for reviews) are a group of seven 
nearby (<lkpc) neutron stars with low (~ 10 32 erg s" 1 ) 
X-ray luminosities and long (3-10 s) spin periods. When 
first discovered, a range of possible reasons for their 
long periods was suggested, with a corresponding range 
in magnetic field strengths of 10 10 — 10 1 5 G : accre- 
tors (iKonenkov fc Popovlll997tlWandll997l). long-period 
pulsa r s" (iKulkarni fc van Kerkwijkl 119981 : iKaplan et alJ 
120021: IZane et alJ l2002h . and middle-aged magnetars 
(|Hevl fe Kulkarnil [l998). Field measurements can allow 
one to distinguish between these possibilities, and thus 
help understand the INS and their place in the greater 
neutron-star population. The magnetic field is also a nec- 
ess ary component of any realistic thermal emission mod - 
els (jZane et alJl2004t iMotch et al.ll2003t iHo et al.ll2007h . 
required to interpret the surface emission and deduce 
radii and other properties of the INS. 

Dipolar magnetic field strengths can be estimated from 
coherent timing solutions, and we used X-ray observa- 
tions to derive such solutions for three INS, 3 finding mag- 
netic fields of (1 — 3) x 10 13 G ( Kaplan fc van Kerkwijkl 
2005a, b, hereafter KvK05a,b; Ivan K erkwiik & Kaplan! 
2008, hereafter lvKK08t also see Ivan Kerkwijk et alJ 
2007|) . This breakthrough is complemented by the dis- 
covery of broad absorption features at energies of 0.2- 
0.75 keV in the thermal (with bolometric luminosities 
of ~ 100 times the spin-down luminosity, the emission 
is certainly thermal) spectra of six of the seven INS. 
Assuming a pure hydrogen atmosphere, fields of 10 - 
10 14 G, and temperatures around 10 6 K, the absorption 
could be due to either proton cyclotron or transitions 
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October. 



between bound states of neutral hydrogen. Intriguingly, 
for the two objects with both spectroscopic and spin- 
down magnetic fields, the agree ment appeared to be good 
(|vaii Kerkwiik fc Kapla"r]|2007l ). 

Here, we constrain the spin-down rate and hence 
magnetic field strength of the INS RX J2143.0+0654 
(also 1RXS J214303.7+065419 or RBS 1774; hereafter 
RX J2143) with a series of dedicated XMM-Newton ob- 
servations. RX J2143 was i dentifi ed as a possible neu- 
tron star by IZampieri et al.l (|200lD on the basis of a soft 
thermal spectrum an d the absence of an optical coun- 
terpart. Using XMM, IZane et al.l (|2005l hereafter IZ+05T ) 
confirmed that the spectrum was soft and blackbody- 
like. They also found a broad absorption feature around 
0.75 keV and identified a candidate 9.437 s periodicity. 
For hydrogen atmospheres, the high energy of the absorp- 
tion (conside rabl y higher than those of the other INS; see 
IHaberj 120071 and Ivan Kerkwiik fc Kaplanl l2007f ) should 
translate to a high magnetic field of > 10 14 G, which 
makes this object a particularly good test for models of 
the origin of the X-ray absorption and for understanding 
the INS population. 

2. OBSERVATIONS AND ANALYSIS 

We observed R X J2143 11 times with XMM 
(|Jansen et al.ll2001f ) in 2007 and 2008, and focus here on 
the data taken with the European Photon Imaging Cam- 
era (EPIC) with pn and MOS detectors, all used in small 
window mode with thin filters (Table[T]). All our observa- 
tions, as well as the one from lZ+05l (taken with the same 
settings), were processed with SAS version 8.0. We used 
epchain and emchain and selected source events from a 
circular region of 37'.'5 radius with energies below 2 keV 
(where flares are negligible; the source is not detected 
above 2 keV) . As recommended, we included only one- 
and two-pixel (single and double patterns 0-4) events 
with no warning flags for pn, and single, double, and 
triple events (patterns 0-12) with the default flag mask 
for MOS1/2. We barycentered the event ti mes using 
the C handra X-ray Observatory position from Rc a et al.l 
(l200l : a = 21 h 43 m 03?38 and 8 = +06°54'17'.'5 (J2000). 
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TABLE 1 

Log of Observations and Times of Arrival 







T7 . . , . a 

hixp. 




1 UA 


Kcv. 


Date 


(ks) 


Counts 


I^IVIJD ) 


820 


2004 May 31 


30.0 


50,032 


53,156.821900(2) 


1360 


2007 May 13 


12.9 


22,915 


54,233.762142(4) 


1362 


2007 May 17 


13.0 


22,926 


54,237.950712(4) 


1368 


2007 May 30 


7.0 


12,369 


54,250.632660(4) 


1375 


2007 Jun 12 


8.3 


14,497 


54,263.916974(5) 


1447 


2007 Nov 03 


8.5 


15,049 


54,407.453138(4) 


1449 


2007 Nov 07 


11.1 


19,480 


54,411.246902(4) 


1449 


2007 Nov 08 


9.5 


16,687 


54,412.213081(3) 


1457 


2007 Nov 23 


5.3 


9,229 


54,427.463785(6) 


1465 


2007 Doc 10 


7.7 


13,517 


54,444.074930(4) 


1538 


2008 May 03 


11.5 


20,812 


54,589.198404(4) 


1546 


2008 May 19 


7.7 


13,414 


54,605.222248(5) 



Note. — All observations used EPIC-pn in the small win- 
dow mode and with the thin filter for both EPIC-pn and EPIC- 
MOS1/2. 

a The exposure time and number of counts given here are for 
EPIC-pn only. b The TOA is defined as the time of maximum 
light of the pulsation closest to the middle of each observation 
computed from the combined EPIC-pn and EPIC-MOS1/2 
datasets, and is given with 1-cr uncertainties. 

2.1. Timing Analysis 

We first determined the energy range that maximized 
the power in a Z\ periodogram for the combined EPIC 
data from 2004 (the longest single observation). We 
found that events below 130 eV (~ 23% of all events) 
were only marginally pulsed, and that removing those 
gave maximum power (Z\ = 58.2; all events gives 
56.9). With these, we found a best-fit frequency of 
v = 0.106066 ± 0.000003 Hz (our final timing solution 
does not change significantly when we include the low- 
energy events). 

Using the above frequency, we determined the times 
of arrival (TO As; see Table [T]) for the combined EPIC 
data from each observation by f itting sin e functions to 
the binned lightcurves (following KvKOSa). The spacing 
and precision of the TOAs is insu fficient fo r an unam- 
biguous timing solution (unlike in lKvK05al b but as in 
vKKOSj). Instead we searched for possible coherent tim- 
ing solutions by iteratively trying sets of cycle counts 
between the different TOAs from 2007 and 2008 (similar 
to IvKKM although we did not incorporate frequency 
information here because it does not add extra informa- 
tion). We find one solution (Fig. [T|) that is considerably 
better than the alternatives (Fig. [2]) with x 2 — 6.0 for 8 
degrees of freedom (dof) and a small, marginally signif- 
icant spin-down rate of v = (—4.6 ± 2.0) x 10~ 16 Hzs -1 
for a dipole magnetic field of J3di P ~ 2 x 10 13 G (see Ta- 
ble[2]). The next best solution has \ 2 = 18.7 and a much 
larger spin-down rate: (—4.7 ± 0.2) x 10 _15 Hzs _1 , re- 
sulting from a difference of 1 cycle for the gap between 
Revs. 1465 and 1538. For a fit with 3 free parameters, 
a change in \ 2 of 12.7 means that the best-fit solution 
is favored at 99.5% confidence. Even less likely solutions 
are found for combinations of other cycle counts (Fig.[2|). 

We confirmed this solution in two ways. First, we ver- 
ified that an analysis using the TOAs from the individ- 
ual instruments gave the same result. Second, we did 
a coherent Z\ periodogram as a function of both v and 
v. The best four solutions identified were the same as 
those found in the TOA analysis. The best-fit peak has 
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Fig. 1. — Phase residuals for RX J2143. In the top panel, 
the residuals are shown relative to a linear model (z> = 0). The 
solid line and filled squares show the best-fit solution including the 
2004 data, while the dashed lines and open squares show the five 
alternate solutions that are within 2<r and have cycle counts for 
the 2004 observation varying by 1 cycle (see text and Fig. [2j . In 
the bottom panel, the residuals to the best-fit quadratic model are 
shown. 

Z\ = 202.9, consistent with the pulsed fraction of ~4%. 
Spin-down is only marginally detected at about 93% con- 
fidence, consistent with the TOA analysis. 

Unfortunately, we cannot unambiguously extend the 
solution back to the 2004 observation: given its offset 
At = 1, 233 days from the reference time to of the above 
solution, the uncertainty on the cycle count is approx- 
imately i<7,>Ai 2 = 1.21 cycles. As a result, about six 
solutions are within ±2cr of the best-fit solution, with 
implied spin-down values differing by 1.21 cycles = 
1.7 x lO^Hzs- 1 (Figs. [Hand©. 

2.2. Spectral Analysis 

We examined the EPIC-pn spectra of RX J 2143 fr om 
our new data to verify whether the basic fits of IZ+05I are 
still valid with our ~ 2.5 times longer total exposure time, 
and to look for po s sible long-term variability such as 
what lde Vries et all (l200l found for RX J0720.4-3125. 
(A full spectral analysis, including the EPIC-MOS and 
RGS data and phase-resolved fits, is in progress.) For 
our spectra, we used the same extraction regions as for 
the timing analysis, and an offset circular region with 
the same radius for the background (which is low). We 
created response files, and binned the spectral files such 
that the bin width was at least 25 eV (about one third 
of the spectral resolution) and the number of counts was 
at least 25. 

We first fit our 11 new observations with an ab- 
sorbed blackbody model over the 0.2-1.5 keV range 
(we use sherpa a nd the xstbabs absorption model of 
[Wilms et al.ll200Cj| ). The best-fit model had column den- 
sity N H = (2.28 ± 0.09) x 10 20 cm- 2 , effective temper- 
ature fcTgfj = 104.0 ± 0.4 eV, and blackbody radius 
i?§B = 3.10 ± 0.04(d/500pc)km, where the blackbody 
parameters are those measured by a n obser ver at infin- 
ity. The results are similar to those of lZ+05l for the 2004 
data, with the slight offset likely resulting from changes 
in the EPIC calibration and from the inclusion of the 
MOS data (for the reprocessed 2004 EPIC-pn data, we 
find values much closer to those given above). The fit 
is reasonable, with \ 2 = 469.0 for 431 dof. The resid- 
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TABLE 2 
Measured and Derived Timing 
Parameters for RX J2143. 0+0654 
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Fig. 2. — Possible timing solutions for RX J2143. Shown are x 2 
values vs. spin-down rates \u\ (along the bottom axis) and dipole 
magnetic field (top axis) with 1 a uncertainties for solutions with 
v < (circles) and v > (squares), using only the data from 
2007-2008 (8 dof). Also shown arc the \ 2 values for the v < 
solutions incorporating the data from 2004 (9 dof; points), which 
have aliases with spacing 1.7 X 10 — 16 Hzs _1 around each of the 
primary solutions (as labeled). 

uals, however, appear to show some systematic devia- 
tions, with dips near 0.75 keV and 0.4 keV. Overall, the 
absorbed and unabsorbed fluxes in the 0.2-2.0 keV band 
are 2.7 x 10~ 12 and 4.0 x 10~ 12 ergs -1 cm" 2 . 

The reduced x 2 for the individual observations ranged 
from 0.9 to 1.5, with w 35 dof each, suggesting the spec- 
tra are consistent with each other. Indeed, fitting for 
both the blackbody temperature and normalization of 
each observation (including that from 2004), but keeping 
ATh fixed, does not decrease x 2 significantly: Ax 2 = 13 
for 22 fewer dof. The changes in temperature found 
are between —0.9% to +1.3%, comparable to the typi- 
cal 2% uncertainties for the individual observations. The 
changes are anti-correlated with those in normalization, 
with i?g B oc kT^, i.e., conserving flux. Thus, we 
see no evidence for variability, and can confidently ex- 
clude the types of changes observed for RX J0720.4— 3125 
(|de Vries et al.ll200l iVink et alJl200l . 

Returning to fits to the 2007-2008 data with a sin- 
gle model, we examined whether the broad absorption 
feature at 0.75 keV was present in our new observations. 
We modeled the absorption by a multiplicative Gaus- 
sian. First, constraining the absorption to be centered at 
0.75 keV w ith a fu ll width at half maximum of 0.06 keV, 
as used in IZ+05L we find a best-fit fractional depth of 
^4 a bs = 22% ±3%, with the remaining parameters similar 
to the blackbody fit and x 2 = 408.2 for 430 dof. Start- 
ing with the same values, but now fitting for the central 
energy and width, we find essentially identical results: a 
best-fit energy of 0.756 keV and width of 0.073 keV, and 
X 2 = 407.0 for 428 dof. A reduction of 61 in x 2 is highly 
significant: an F test 4 gives a probability of random oc- 
currence of 4 x 10~ 13 . To verify this, we simulated 1000 
sets of 11 observations like ours using the best-fit black- 
body parameters, and fitted those to blackbody models 
with and without absorption. We found that the distri- 
bution of j4 a bs was centered at zero with a root-mean- 

4 Using an F test for finding lines is incorrect wh en the line model 
is loc ated at the boundary of parameter space (Protassov ct al. 
2002?). This, however, is not the case here since A a b s can have 
either sign. 



Quantity 



Value 



Dates (MJD) . . 54,234-54,605 

t (MJD) 54,383.648930(2) 

v (Hz) 0.1060644595(11) 

v (lO" 16 Hzs" 1 ) -4.6(20) 

TOA rms (s) ... 0.3 

X 2 /dof 6.0/8 

P (s) 9.42822888(9) 

P (lO-^ss" 1 ). 4.1(18) 

Tchar (Myr) 3.7 

Bdip (10 13 G) . . 2.0 

E (10 30 ergs" 1 ) 1.9 

Note. — Quantities in parentheses 
are the formal 1-cr uncertainties on the 
last digit. T c h ar = P/2P is the char- 
acteristic age, assuming an initial spin 
period Pq -C P and a constant mag- 
netic field; B dip = 3.2 X W w sJ~Pp 
is the magnetic field inferred assum- 
ing spin-down by dipole radiation; E = 
3.9 X 10 46 P/P 3 is the spin-down lumi- 
nosity. 

square spread of 3%, similar to the uncertainty on A a i, s 
in our fit to the actual data. The distribution of x 2 also 
conforms to expectations, with the x 2 including absorp- 
tion never differing from the blackbody x 2 by more than 
8.5 (compared to 61 for the real data). Th us, we c onfirm 
the detection of absorption at 0.75 keV bv lZ+05l Using 
the same continuum model but letting the absorption 
depth vary for each observation, we again find no sta- 
tistically significant variability (Ax 2 = 10 for 11 fewer 
dof). 

If we vary the energy of the absorption line -E a b s over 
the range 0.3-1.0keV (avoiding the edges of the spec- 
trum, where the continuum fit and the absorption are 
highly covariant), we find that, with one exception, only 
for i? a bs — 0.75 keV is there significant absorption. Oth- 
erwise, Ax 2 < 10 and A a ^ s is around 0. The one ex- 
ception is a hint of absorption at 0.42 keV, which gives 
Ax 2 = 15 and ^ abs = 9.7% ± 2.3%. This detection is 
significant at the ~ 3er level given the 15 trials that we 
did (none of our 1000 simulations achieved such a high 
Ax 2 )- Such a line also improves the fit for the 2004 
data, although not by a statistically significant amount 
(Ax 2 = 4). Confirmation of this with the EPIC-MOS 
and RGS data is ongoing. 

3. DISCUSSION AND CONCLUSIONS 

A large magnetic field of > 10 14 G is inferred 5 for 
RX J2143 if one interprets the 0.75 keV absorption in 
its spectrum as arising from either proton cyclotron 
absorption (1.4 x 10 14 G for a gravitational redshift 
z = .3) or io nization of hydrogen (even hi gher, as in 
IZ+05t also see Ivan Kerkwiik fc Kaplanll2007l ). This sce- 
nario appeared to work well for RX J0720.4— 3125 and 
RX J1308. 6+2127, where those transitions could match 
the observed absorption features for magnetic fields of a 

5 This interpretation ignores the suppression of absorption lines 
due to va cuum resonance m ode conversion in such strong magnetic 
fields; sec IHo fc Lail J2004T) . 



4 



Kaplan & van Kerkwijk 



o 

w 3.5 

T) 

E 3 

U 

1 2.5 



j 




1 1 
1 1 

/ " 


- 


I 




1 1 

1 1 
1 1 
1 1 
1 1 


1308 


1 




1 1 
1 1 

1 1 
1 1 


• 

□ 


1 / 




0720 [[ 

• .. 

1 1 
1 1 


□ 

if 143 


y 


1856 
• 


1 1 
1 1 
1 1 
1 1 
1 1 
1 1 


□ 


o[ 
o 1 




1 1 

5 1 o ~ 





1.5 1 |J 1 1 " 1 1 

40 60 80 100 120 140 160 



Effective Temperature (eV) 

Fig. 3. — Dipole magnetic field vs. effective temperature for the 
seven INS. Temperatures are derived from blackbody fits including 
absorption features, and have been correcte d for a nomina l grav i- 
tational redshi f t of 0. 3. They are taken from Burwitz ct al. (2003), 
Ide Vries et alJ (1200411 . Ivan Kerkwijk et al.l ((2004 ). ISchwope et at] 
(200 5), Z+Q5| (which is c o nsiste nt with the results from Sec- 
tion |5j5), and Habcrl ct al. (2004). The magnetic fields are those 
estimat ed from spin-down, and are taken from IKvK05al .b and 
vKK08|; sources without spin-down solutions are shown as ver- 
tical lines. For RX J2143, the main solution from Table [2] is 
shown, as well as the possible aliases including the 2004 data 
Fig-O • The diagonal line is where a surface of iron would condense 
Medin fc Lai 200?j) . Note the importance of using the spin-down 
derived magnetic fields; the inferences made in less secure ways can 
be drastically diff erent (compare to Fig. 1 in Turolla ct al. 2004 and 
IPons et al.ll2007D . 

few tim es 10 13 G, comparable to what was inferred from 
timing (|van Kerkwijk fc Kaplanll2007t lHaberlll2007l and 
references therein). It also seemed consistent with the 
lack of features in RX J1856.5— 3754, since that source 
has the weakest field (|vKK08t ). 

For RX J2143, though, the model breaks down: the 
required strong magnetic field is inconsistent with our 
timing measurements (-Bdip = 2 x 10 13 G) at the 10~ 4 
level (given A\ 2 > 20 for 3 parameters). In principle, 
the possible absorption feature at 0.4 keV may be eas- 
ier to accommodate, although it still occurs at higher 
energy than the features seen in RX J0720.4— 3125 and 
RX J1308. 6+2127 while the magnetic field that we mea- 
sure here is nominally weaker than the fields of those 
sources (Fig. [3]). Furthermore, this leaves the 0.75 keV 
feature unexplained. It being a "harmonic" of the 0.4 keV 
line seems unlikely, as it is almost twice as strong, 
while o ne would expect harm onics to be significantly 
weaker (jPavlov fc Panovi 119761 : G. G. Pavlov 2007, pri- 
vate communication). A better match may be possible 
with a n atmosphere with helium or even heavier ele- 
ments (iPavlov fc Bezchastnovl2005l : lHailev fc Morill2002t 
IZ+Ol iMori fc Hdl2007l) : of course, in this case hydrogen 
must be absent, given the rap id gravitational settling 
time (|Alcock fc Illarionovll 19801) and the small amount of 
hydrogen required to be optically thick (|Romanilll987t ). 

Geometry may offer a partial solution. The magnetic 
field that we infer from the spin-down rate is actually 
the true dipolar field strength times a function of the 
angle a between the magnetic and rotation axes. That 
function is si n a in the traditional vacuum dipole model 
(|Pacini|[l967l ). which would give a large range of possi- 
ble true field strengths. More modern analyses, however, 

give something more like vT + sin 2 a (ISpitkovskvll2008l) . 
and thus a range of only 40%. Therefore, it seems un- 



likely that this can remedy the situation. Modeling of 
the pulse profile and phase-resolved spec troscopy has the 
poten tial to constrain the geometry (e.g. JZane fc TuroTIal 
2006), although X -ray polarimetry may be required for 
unambiguous results. 

Turning to the properties of the INS as a whole, 
our result shows that not only do the INS clus- 
ter at long periods, but they are close in both 
axes of the P — P plane. The inferred mag- 
netic fields for RX J1856. 5-3754, RX J0720.4-3125, 
RX J1308. 6+2127, and now RX J2143 are all within a 
factor of 2 of each other. Intriguingly, we also find a pos- 
sible correlation between the effective temperature and 
m agnetic field (Fig . [3j a similar realization was made 
bv lPons et al.ll2007| for a larger but less uniform sample, 
but see also lTurolla et al.ll2004j where the use of magnetic 
fields inferred from spectroscopy led to different conclu- 
sions). Whether or not RX J2143 fits this trend, however, 
will require an unambiguous timing solution. 

A possible explanation for the clustering of the mag- 
netic field strengths, and perhaps also the correlation 
with temperature, is that the fields were originally sig- 
nificantly stronger (10 14 -10 15 G), and decayed. Such a 
scenario was originally proposed as a way to keep the INS 
hotter for muc h longer and thus ma ke the n atal popula- 
tion s maller (jHevl fc Kulka rnlll998t also see lColpi et alJ 
2000). It is unlikely that field decay contributes much 
to the current thermal state, as the cooling and kine- 
matic ages (from tracing the INS back to birth locations) 
agree well and the products of those ages and X-ray lu- 
minosities are about 100 times the energy currently in 
the dipole magnetic field, but decay may still have influ- 
enced the current field strengths. In particular, relatively 
fast field decay naturally leads to the tightly grouped 
perio ds and magnetic fields of the INS (jPons fc GeppertJ 
|2007| ). It also naturally resolves why the spin-down ages 
P/2P (which assume constant magnetic field) are sig- 
nificantly lo nger than the kinematic ages: ~2Myr ver- 
sus . 6 My r (Kaplan ct al. 2 0071 : Ivan Kerkwijk fc Kaplanl 
120071 : lvKK08h Finally, field decay might also lead to 
a correlation between field and temperature, since both 
would be a function of age (note, however, that based 
on kinematic ages, RX J1856.5— 3754 appears younger 
than RX J0720.4— 3 125 even though the former is colder; 
Kaplan et ail 120071 ). If field decay occurred, the INS 
would be the descendants of something like the magne- 
tars, having had stronger fields (>2 x 10 14 G, especially 
in the interiors) in the past, but with merely above aver- 
age fields now. It could point to an evolutionary differ- 
ence between the INS and the h igh-magnetic- field radio 
pulsars (e.g.. ICamilo et aLll2000D that inhabit the same 
region of the P-P diagram. 

The correlation between temperature and field 
strength might also be evidence for something rather 
different, viz., surface physics. For a given magnetic 
field, as the surface cools, it is expected to condense 
(jRudermanl Il97i lLall r200lh . For iron, the condensa- 
tion line is quite similar to the correla tion we observe 
(|Turolla et alJ [200l iMedin fc La3 [2007h . From Fig. p 
one sees there is about a factor 2 difference in magnetic 
field, but this may just reflect the difference between the 
dipole and crustal fields, or between the true temper- 
ature and that inferred from a blackbody fit (indeed, 



The Spin-down of the Nearby Isolated Neutron Star RX J2143.0+0654 



5 



differences in the r ight sense are expected from atmo- 
sphere models; see IZavlin fc Pavlov! l2002f ) . If conden- 
sation indeed plays a role, it would help determine the 
surface composition (lighter elements do not condense 
as easily), perhaps help understand the peculiar spectra 
of the INS, and, since a conde nsed surface may i nhibit 
the formation of a vacuum gap (jMedin fc Lail l2007l. help 
explain the lack of rad io emission from the INS (e.g., 
iKondratiev et all [2008). It is less clear, however, what 
would keep the INS on the condensation line, or make 
them evolve along it. 

The above possibilities may well help understand 
what separates the INS from "normal" rotation-powered 
middle-aged pulsars, which have a much wider range of 
magnet i c field strengths in a comparable sample (e.g., 
iKaplanl I2008T) . Understanding how the fields of these 
sources evolve and how this couples to the surface tem- 



perature is also required to de rive meaningful c onstraints 
from cooling measurements ()Page et al.l I2004D . To see 
whether the fields are indeed clustered and/or correlated 
with temperature, will require measurements of other 
sources, and refinement of that of RX J2143. 



Based on observations obtained with XMM-Newton, 
an ESA science mission with instruments and contribu- 
tions directly funded by ESA Member States and NASA. 
DLK was supported by NASA through Hubble Fellow- 
ship grant #01207. 01-A awarded by the Space Telescope 
Science Institute, which is operated by the Association of 
Universities for Research in Astronomy, Inc., for NASA, 
under contract NAS 5-26555. 
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